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Abstract 

Stop particles are expected to be the lightest squarks in supersymmetric theories and 
the search for these particles is an important experimental task. We therefore present 
the cross sections for the production processes pp/pp — ► t±ti and £2^2 at Tevatron 
and LHC energies in next-to-leading order supersymmetric QCD. The corrections 
stabilize the theoretical predictions for the cross sections, and they are positive, 
thus raising the cross sections to values above the leading-order predictions. Mixed 
tit2/t\t2 pairs can only be generated in higher orders at strongly suppressed rates. 
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1 Introduction 



Within the squark sector of supersymmetric theories, the top-squark (stop) eigenstate t\ is 
expected to be the lightest particle Q. If the scalar masses in grand unified theories, for 
instance, are evolved from universal values at the GUT scale down to low scales the t\ top- 
squark drops to the lowest value in the squark spectrum. Moreover, the strong Yukawa coupling 
between top/stop and Higgs fields gives rise to large mixing, leading to a potentially small mass 
eigenvalue for i\ @. 

In e + e~ and pp/pp collisions stop particles are produced in pairs. Present limits from LEP2 
indicate a t\ mass in excess of 67 GeV, independent of the mixing angle in the stop sector, but 
depending on the lightest neutralino mass Q . Preliminary analyses at the Tevatron have led to 
a lower limit of 93 GeV for the t\ mass, depending on the lightest neutralino mass @. 

The cross sections for the production of squarks (q) and gluinos (g) in hadron collisions 
have been calculated at the Born level already quite some time ago ||. Only recently have 
these theoretical predictions been improved by calculations of the next-to-leading order (NLO) 
SUSY-QCD corrections for squark/gluino production, with the final-state squarks restricted to 
the light-flavor sector (g ^ t) |J. In the present paper we supplement this analysis by the 
corresponding analysis for the stop sector 

pp/pp — ► titi + X and t 2 t2 + X (1) 

in the pp collisions of the Tevatron and the pp collisions of the LHC. This NLO calculation 
is motivated by two requirements: First, to stabilize the theoretical predictions for the cross 
sections with respect to the renormalization and factorization scales, which introduce spurious 
parameters into fixed-order calculations. And second, to improve the accuracy of the theoretical 
predictions for the cross sections. Extending the light-flavor analysis to stop production is a 
necessary step since the mixing effects must be taken into account properly. At NLO [i.e. 0(aP s )] 
the production cross sections are still diagonal in the stop sector. The production of mixed 
t-\tz/t$b\ pairs is suppressed as the cross section is of order af. Since the calculation is rather 
involved but the rate is small, we have studied this process in the limit of large gluino mass to 
exemplify the expected size of the non-diagonal cross section. 

The cross sections for the diagonal production of stop particles, Eq. (||), depend essentially 
only on the mass of the produced stop particles, raj or . The dependence of the cross sections 
on the other SUSY parameters, i.e. the gluino mass nig, the masses of the other squarks and 
the mixing angle 0, is very weak since these parameters affect only the higher-order corrections 
and are not relevant at leading order. 

The results for the diagonal production of stop particles, t\t\ and £2^2 ; will be discussed in 
the next section. The size of the cross section for mixed-pair production txt^ and will be 
estimated subsequently. 
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2 Diagonal stop-pair production 



At hadron colliders, diagonal pairs of stop particles can be produced at lowest order QCD in 
quark-antiquark annihilation and gluon-gluon fusion: 



qq — * t±t± and £2*2 
99 ~^ t\t\ and £2^2 



(2) 



Mixed pairs t-^t^ an d t<tf,\ cannot be produced in lowest order since the gtt and ggtt vertices are 
diagonal in the chiral as well as in the mass basis. The relevant diagrams for the reactions @ 
are shown in Fig.|l|a. The corresponding cross sections for these partonic subprocesses may be 
written as (k = 1, 2): 



&Lo[qq -> hh\ 



a 2 a 7r 2 



27 



2 

a Q 7r 



h 48 + ~W 




(3) 
(4) 



The invariant energy of the subprocess is denoted by y/s, the velocity by 0k = wl — 4m? / s. 

The cross sections coincide with the corresponding expressions for light-flavor squarks in the 
limit of large gluino masses, cf. Ref. ||. 

The hadronic pp/pp cross sections are obtained by folding the partonic cross sections with 
the qq and gg luminosities. At the Tevatron the dominant mechanism for large stop masses is 
the valence qq annihilation. The fraction of qq events rises from 0.55 to 0.86, if the t\ mass is 
increased from 100 to 200 GeV. At the LHC the gluon-fusion mechanism plays a more prominent 
role. For a t\ mass below 200 GeV, more than 90% of the events are generated by gg fusion, 
cf. Table @. 

We closely follow the approach of Ref. || for the calculation of the NLO SUSY-QCD correc- 
tions. The virtual O (a s ) corrections involve the usual SUSY-QCD corrections to the propaga- 
tors and vertices, as well as box and rescattering diagrams, cf. Ref. ||. To this order the mixing 
angle 9 enters the cross section only through corrections involving the ttg and the four-squark 
couplings. The relevant couplings of this type are described by the Lagrangeans 



£3 



V2g s TtAl + V 12 )h 



cos 9 hi -75) - sin 6*^(1 +75) tjii* + h.c. 



(1 + P12) hi hj { cos 2 (2#~) Af l h* k in + 2 sin 2 (20) A l f l - A l ( 



iki 



jki 



with the color tensors (iV c =3) 



kl 



A\ j 



j^jkl 



+ 4 cos(20) Af l £ (QLt QLi ~ qui Qm) 



5 il 5 kj -5 ij 5 kl /N c 
5 il 5 kj + 6 ij 5 kl ) (N c - 1)/N C 



(5) 



(6) 
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Here q^lR represent the left(L)- and right(R)-chirality light-flavor squarks. Indices of the fun- 
damental/adjoint representation of color SU(3) are denoted by i, j, k and a, respectively, and 
the generators of the fundamental representation by T a . The operator V12 permutes the 1- and 
2-components of the top-squarks: 

V12 '■ [ti t2) cos 9 — > — sin9, sinfl^cosc?] (7) 

Two typical diagams in which 9 enters the scattering amplitude gg — > tktk through vertex and 
rescattering corrections are shown in Fig. [l]b. Since mixing enters explicitly only through higher- 
order diagrams, the angle 9 need not be renormalized [Q] in the present calculation and it can 
be identified with the lowest-order expression derived from the stop mass matrix. 

As usual, the virtual corrections are supplemented by gluon radiation from color lines and 
vertices, as well as contributions from the inelastic Compton processes qg — ► tktkq and gq — > 
tktkq, cf. Ref.jg] for diagrammatic details. 



The singularities associated with the NLO corrections are isolated by means of dimensional 
regularization and renormalized within the MS scheme.Q The renormalization of the QCD 
coupling is performed in such a way that the heavy particles (top quarks, gluinos, and light- 
flavor squarks) are decoupled smoothly for momenta smaller than their masses. This implies that 
the heavy particles do not contribute to the evolution of the couplings and parton densities. The 
masses of the light quarks (q 7^ t) are neglected and the top-quark mass is set to m t = 175 GeV. 
For the calculation of gluon bremsstrahlung, the phase space for gluon radiation is split into 
two distinct regimes, one accounting for soft gluons and the other for the hard gluons. The 
separation is implemented by introducing a cut-off parameter A in the invariant mass of the 
radiated gluon and one of the heavy particles in the final state. The cut-off parameter is chosen 
so small that it can be neglected with respect to any other mass scale in the process. As a result 
of the split-up of the phase space, terms of the form log* A (i = 1,2) occur in both the soft 
and the hard cross sections. If soft and hard contributions are added up, any A dependence 
disappears from the cross sections in the limit A — > 0. 

At lowest order, the cross sections for diagonal t\t\ and ^2 pair production are given by the 
same analytical expression. At next-to-leading order, the ttg and four-squark interactions will 
affect the production cross sections for the two diagonal pairs in different ways, introducing the 
explicit dependence on the mixing angle. However, the 9 dependence will turn out to be very 
mild. It follows from Eqs. (5), that the analysis of £2*2 pairs can be copied from the analysis of 
t\t\ pairs after applying the permutation V12 defined in Eq. (fjj). 

For the detailed description of the partonic cross-sections we introduce scaling functions /, 



ra- 
ti 



fZ +S {ri,r,~9) + f"iv) + hmo. 

\ in- , 




(8) 

with r generically denoting all possible mass ratios m/mf . The indices i,j = g,q,q indicate 
the partonic initial state of the reaction ij — ► if-tp.- The center-of-mass energy of the partonic 



1 Note that the spurious breaking of supersymmetry in the MS scheme |j] has no effect on the NLO stop-pair 
cross sections. The Yukawa couplings become effective only as a part of the higher-order corrections and need 
not be renormalized. 
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reaction, y/s, is absorbed in the quantity r\ = s/4m~ — 1, which is better suited for analyzing 
the scaling functions in the various regions of interest. We have identified the renormalization 
and factorization scales: hr = \ip = /i. The scaling functions are divided into the Born term 
f B , the sum of virtual and soft-gluon corrections f v+s , the hard-gluon corrections f H , and 
the scale-dependent contributions /. The log* A (i = 1,2) terms are separated from the soft- 
gluon corrections and added to the hard-gluon part. The hard-gluon corrections are therefore 
independent of the cut-off for A —* 0. 

The scaling functions are displayed in Fig.|| for the quark-antiquark, gluon-gluon, and 
(anti)quark-gluon channels. The mixing angle and the gluino and light-flavor squark masses are 
defined in a minimal supergravity scenario that will be discussed later in detail. The dependence 
on these parameters is weak. In particular if the mixing angle is varied over the full range, the 
impact on f v+s , the only scaling function affected by 9, is very small. Due to the P 3 behavior 
of the LO quark-antiquark cross section, not much structure is observed in this channel in 
the vicinity of the production threshold (/3,rj -C 1). However, for the gluon-gluon channel two 
sources of large corrections can be identified in the threshold region. First, the exchange of (long- 
range) Coulomb gluons between the slowly moving massive particles in the final state leads to 
a singular Sommerfeld correction ~ ira s / f3, which compensates the LO phase-space suppression 
j3. The scaling function f^ g +s therefore approaches a non-zero value at the threshold. [For 
the qq incoming state the cross section remains suppressed ~ /3 2 , since the LO gluino-exchange 
contribution that dominates the threshold behavior of light-flavor squark-pair production |J, 
is absent for stop-pair production]. It should be noted, however, that the screening due to the 
non-zero lifetimes of the top-squarks reduces the Coulomb effect considerably. Second, as a 
result of the strong energy dependence of the gluon-gluon cross sections near threshold, large 
gluonic initial-state corrections of the type /31og*/3 (i = 1,2) emerge. The leading log 2 /5 terms 
are universal and can be exponentiated. Near threshold, the scaling functions can be expanded 
in /?: 



J 99 384 Jqq 54 

fV+S _ fB 11 f V_+S _ _fB_ 1 

J 99 J 99 3mf3 Jqq J ^48/3 



fH _ fB 

■>gg hg 



fH _ fB_ 
Jqq Jqq 



^log 2 (8/3 2 )-^log(8/3 2 



hg = -fgg MS/? 2 ) f qq - = -f* l g(8/? 2 ) 



2 

2vr 2 ~ m ^ ' Jqq Jqq 3^ 2 

At high energies the NLO partonic cross sections in the gluon-gluon and (anti)quark-gluon 
channels approach non-zero limits asymptotically, to be contrasted with the scaling behavior 
~ 1/s of the LO cross sections. This is caused by the nearly on-shell exchange of space-like 
gluons, associated with the inelastic Compton processes and gluon radiation in the fusion process. 
Exploiting the factorization in transverse gluon momentum at high energies, the high-energy 
scaling functions can be determined analytically: 

f H = f H = ^_ (1Q) 

J " 43200^ Jqg 194400tt 1 ' 

f =-— f — 
J " 720tt Jqg 3240vr 
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The scaling functions f gq - are identical to f qg . The ratio of the f gg and f qg scaling functions is 
given by 9 : 2, corresponding to the probability for emitting a soft gluon from a gluon (twice) 
or a quark. 

The numerical analyses of the hadronic cross sections have been performed for the Fermilab 
Tevatron pp collider with a center-of-mass energy of y/S = 1.8 TeV, and for the CERN Large 
Hadron Collider (LHC) with a pp center-of-mass energy of yfS = 14 TeV. We have adopted 
the CTEQ4M parametrization of the parton densities Q. The uncertainty due to different 
parametrizations of the parton densities in NLO is less than ~ 5% at the Tevatron and less than 
~ 10% at the LHC (irrespective of the choice for the scale //). The difference between the two 
estimates can be attributed to the fact that the experimentally well-determined valence quarks 
dominate at the Tevatron, while the gluons are more prominent at the LHC. 

In Fig.^ we present the dependence of the total cross sections for titi production on the 
renormalization and factorization scale /i = [ir = f^F- For a consistent comparison of the 
LO and NLO results, we have calculated all quantities [a<j(/^), the parton densities, and the 
partonic cross sections] in LO and NLO, respectively. In LO the scale dependence is steep and 
monotonic: changing the scale from /i = 2m^ to fj, = the LO cross section increases 

by more than 100%. In NLO the scale dependence is strongly reduced, to about 30% in this 
interval at the Tevatron. At the same time the cross section is considerably enhanced at the 
central scale (fi = ). The results are qualitatively similar for the LHC. 

The magnitude of the SUSY-QCD corrections is illustrated by the K factors in Table [l] for 
Tevatron and LHC energies. The K factor is defined as K = (JnloI&lOi with all quantities 
calculated consistently in lowest and in next-to- leading order. In these calculations the scale is 
fixed at the central value (fi = m t k )- For illustration, the super symmetric parameters have been 
fixed within the minimal superg ravity (SUGRA) modelQ such that m g = 284 (627) GeV and 
sin(2(9) = —0.99 (—0.94) for the Tevatron and the LHC, respectively. The stop masses in these 
scenarios are given by m? = 153 (325) GeV and = 347 (592) GeV, all other squarks are 
assumed to be mass degenerate with m g = 256 (539) GeV. However, in order to focus on the 
mass dependence of the cross sections and of the NLO corrections, the mass of the produced stop 
particles is varied around the SUGRA-inspired central value, independently of the other SUSY 
parameters. In the mass range considered, the SUSY-QCD corrections are small (and negative) 
if the qq initial state dominates, cf. Table |l[ If, in contrast, the gg initial state dominates, the 
corrections are positive and reach a level of 30-40%. The relatively large mass dependence of 
the K factors for t\t\ production at the Tevatron can therefore be attributed to the fact that 
the gg initial state is important for small mj , whereas the qq initial state dominates for large 

The total cross sections play a crucial role in the experimental analyses. They either serve to 
extract the exclusion limits for the mass parameters from the data, or, in the case of discovery, 
they can be exploited to determine the masses of the stop particles. The total cross sections for 
pp/pp — ► fit 1,^2^2 are given in Figs. f|/]5] for the Tevatron/LHC, respectively. The cross sections 
depend essentially only on the masses of the produced stop particles, and very little on the other 

2 For an approximate solution of the renormalization-group equations in the SUGRA-inspired model, the pro- 
gram SPYTHIA jl(J has been adopted. For the Tevatron (LHC) parameters the input values mo = 100 GeV, 
mx/2 = 100 (250) GeV, Aq = 300 GeV, tan/3 = 1.75, and /j, > have been chosen, generating the mass and 
mixing parameters quoted above. 
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m- t [GeV] 


-^■Tevatron 


ggm 


QQin 


-Klhc 


99m 


QQin 


hh 


70 


1.41 


0.64 


0.36 


1.25 


0.96 


0.04 




110 


1.30 


0.41 


0.59 


1.32 


0.95 


0.05 




150 


1.19 


0.25 


0.75 


1.37 


0.94 


0.06 




190 


1.11 


0.15 


0.85 


1.40 


0.92 


0.08 


hh 


280 


1.01 


0.06 


0.94 


1.44 


0.89 


0.11 




320 


1.00 


0.04 


0.96 


1.45 


0.88 


0.12 




360 


0.98 


0.03 


0.97 


1.46 


0.86 


0.14 




400 


0.95 


0.02 


0.98 


1.48 


0.85 


0.15 



Table 1: K factors for diagonal stop-pair production at the Tevatron and the LHC for a sample of 
stop masses. Scale choice: [x = m^. The SUGRA-inspired parameters adopted in the calculation 
of the higher- order corrections are defined in the text. For completeness also the LO initial-state 
gg and qq fractions are given. 

supersymmetric parameters, i.e. the gluino mass, the masses of the light-flavor squarks and the 
mixing angle. The variation of the cross section with the gluino mass and the mixing angle is 
indicated by the thick NLO curves. With cross sections typically in the range between 1 and 
100 pb, sufficiently large samples of 10 3 to 10 5 stop events can be accumulated at the Tevatron 
for an integrated luminosity of J C = 1 fb~ , provided these particles exist in the mass range 
below 200 GeV. For an integrated luminosity of f C = 100 fb _1 at the LHC a large sample of 
10 5 to 10 7 stop events could be collected for masses in the 200-500 GeV range. 



3 Mixed t\ti and t^tx pairs 

In pp/pp collisions, the mixed final states hh and hh cannot be produced in lowest order. The 
production cross section for non-diagonal stop pairs is therefore of order . This higher-order 
cross section is small but complicated to calculate. We therefore discuss the size of the cross 
section in the limit of large gluino mass. In this limit, the four-squark couplings relevant for 
non-diagonal stop-pair production are described by the Lagrangean 

C' A = & S m(26) (h* t 2j + h* hj) | cos(2fl) Af l (h* k hi - hi hi) 

+ 2 Af l (&* QLl ~ Qd QRl) } (11) 

[where however the second term does not contribute to the cross section for mass degenerate 
left / right light-flavor squarks q] . Only two one- loop diagrams contribute to the scattering ampli- 
tude in this limit (see Fig.|l|c). They involve the production of diagonal t^h pairs in gg collisions, 
subsequently transformed to non-diagonal t\t2 and t2h pairs by rescattering in the final state. 
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The loops can easily be evaluated: 



Poo[gg -» hh + hh] = sin 2 (40) 



37 



13824 2vrs 3 



log (x{) 



(12) 



where the subscript in the cross section o"oo indicates the limit m g — > oo. The coefficient A 1 / 2 
is the usual 2-particle phase-space factor, i.e. A = [s — (mj + ^ 2 ) 2 ][s — — "it 2 ) 2 ]/s 2 , and 
£fc = (A — + 1); the logarithmic singularities are defined properly by the infinitesimal 

shift s — ► s + is in /3^. 

The cross section depends strongly on the mixing angle 9 through the overall factor sin 2 (40). 
Numerical values for the diagonal and non-diagonal pair cross sections are compared in Table ^ 
for the default SUGRA- inspired SUSY parameter sets fl(| adopted already earlier. Note that 
the mixed-pair cross section is given in this table without the mixing factor sin 2 (40). Evidently, 
the values for the cross section for producing mixed stop pairs in the large nig limit are very 
small at the Tevatron as well as at the LHC. 





<7[fb] 


_ limit 
Oqq U qq 


limit 
Vgg O-gg 


Tevatron 


hh 
hh 

t\t2 + ^2^1 


0.64 -10 3 0.64 -10 3 
1.51 1.54 



0.42 -10 3 0.42 -10 3 
0.105 0.108 

1.54 -lO" 4 


LHC 


hh 
hh 

t\t2 + ^2^1 


0.60 -10 3 0.59 -10 3 
37.1 37.1 



6.50 -10 3 6.47 -10 3 
0.23 -10 3 0.24 -10 3 
3.53 -10" 2 



Table 2: Cross sections for diagonal and non-diagonal pair production at the Tevatron and the 
LHC, using the default SUGRA-inspired values for the SUSY parameters. The non-diagonal 
results are given without the mixing factor sin 2 (40). Scale choice: average mass of the produced 
stop particles. The superscript 'limit' denotes the asymptotic value of the cross section for large 
gluino masses. 



4 Summary 

The picture that has emerged from the SUSY-QCD analysis, is quite simple, (i) The cross 
sections for the production of diagonal pairs pp/pp — > Mi, ^2 depend essentially only on the 
masses of the stop particles produced and very little on the other supersymmetric parameters. 
Bounds on the t\ti production cross section can therefore easily be translated into lower bounds 
on the lightest stop mass without reference to other supersymmetric parameters. On the other 
hand, if stop particles were to be discovered, the cross section can be exploited directly to 
determine the two stop masses and m^. (ii) If mixed stop pairs could be discovered in 

pp — ► hh + hti for sufficiently high integrated luminosity, the mixing angle can be derived 



8 



from the magnitude of the cross section which is proportional to sin 2 (40). In contrast to mixed 
stop-pair production via Z exchange in e + e~ annihilation the cross section is suppressed 
by 0{a 2 s ) with respect to diagonal pair production. 

Acknowledgements: We are very grateful to S. Lammel for advise and discussions on the 
search for stop particles at the Tevatron. Collaboration with R. Hdpker is also gratefully ac- 
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Note: The program for calculating the stop cross sections has been included into the package 



PROSPINO, http : / / wwwcn. cern. ch / ^mspir a . 
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Figure 1: Generic Feynman diagrams for the production of pairs of stop particles: (a) Born 
diagrams for quark-antiquark annihilation and gluon fusion; (b) higher-order diagrams for the 
diagonal production including stop mixing (dotted vertices); (c) non-diagonal production in the 
limit of decoupled gluinos (mixing vertices are dotted). 
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-0.002 1 — — — — — — 1 

-3 -2-1 2 3 

10 10 10 1 10 10 10 

Figure 2: The scaling functions for the production oft\i\ pairs as a function of rj = s/4m~ — 1. 
The notation follows Eq. (^). The variation of the scaling function fij +S for all possible values 

of the mixing angle 9 is indicated by the line-thickness of the corresponding curves. The scaling 
functions f g q are identical to f qg . 
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Figure 3: Renormalization/factorization-scale dependence of the total cross sections for ti -pair 
production at the Tevatron and the LHC. The SUSY mass parameters correspond to the central 
values of the SUGRA-inspired scenario described in the text. 
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Figure 4: The total cross sections for the production of pairs of stop particles (t^tk) at the 
Tevatron as a function of the stop masses. The band for the NLO result indicates the uncertainty 
due to the renormalization/ factorization scale. The light-flavor squark masses, the gluino mass 
and the mixing parameter are derived within the SUGRA-inspired scenario defined in the text. 
The line-thickness of the NLO curves represents the simultaneous variation of the gluino mass 
between 200 (284) an d 800 GeV for t\(t-2)-pair production and the variation of sin(20) over its 
full range. 
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Figure 5: The same as Fig.% but for the LHC. The SUSY mass spectrum is described in the 
text. The gluino mass is varied between 400 (600) and 900 GeV for t\{t^)-pair production. 
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